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ABSTRACT 


Air-sea transfers of sensible heat, latent heat, and momentum are 
computed from twenty-five years of middle-latitude and subtropical ocean 
weather ship data in the North Atlantic and North Pacific using the bulk 
aerodynamic method. The results show that monthly-averaged wind speeds, 
temperatures and humidities can be used to estimate the monthly-averaged 
sensible and latent heat fluxes computed from the bulk aerodynamic equa- 
tions to within a relative error of approximately 10%. The estimates of 
monthly-averaged wind stress under the assumption of neutral stability 
are shown to be within approximately 5% of the monthly-averaged non-neutral 
values. 


1. Introduction 

Estimates of the transfers of sensible heat, latent heat, and momen- 
tum between the atmosphere and the oceans may provide important insight 
into the dynamics of interannual and decadal climate fluctuations. Sensi- 
ble and latent fluxes combine with radiative fluxes to determine the major 
part of the net air-sea exchange. Momentum transfers play an important 
role in driving the ocean general circulation and particularly in deter- 
mining the vertical structure of the upper ocean. 

The first systematic evaluation of the sensible and latent heat trans- 
fers over the world ocean was performed by Budyko (1963). He computed 
long-term monthly-mean values of the near-surface wind speed and of the air- 
sea temperature and humidity differences and evaluated the fluxes by the 
well-known bulk aerodynamic formulas with constant transfer coefficients 
(see section 2). 

Budyko* s method depends on the validity of the following two conditions 

1. The covariance of the wind speed with either the air-sea temper- 
ature difference or the air-sea humidity difference must be less 
than the product of their respective mean values (except for cases 
in which the fluxes are small). This condition was found to be 
satisfied in several important climatic regions of the world 
oceans for limited time periods by Malkus (1962), Budyko and 
Gandin (1966), Kraus and Morrison (1966) and Fissel et al. (1977). 

2. The transfer coefficients for sensible heat and evaporation must 
not be strongly dependent on either wind speed or air-sea temper- 
ature and humidity differences for typical conditions over the 
oceans. Budyko (1974) concludes that the stability of the surface 
boundary layer, as parameterized by the bulk Richardson number, 

is almost neutral except in'weak wind conditions. As we will show 
in subsequent sections, recent developments in surface boundary 
layer theory support Budyko 's conclusion. 

The first systematic evaluation of air-sea momentum transfers on a 
global basis was presented by Hellerman (1967). His computations were 
based on the bulk aerodynamic formulas with a drag coefficient for neutral 
surface boundary layer conditions and wind data consisting of monthly- 
averaged wind roses with both wind direction and speed categories. The use 
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of only one mean wind speed in each direction category of a monthly- 
averaged wind rose leads (see Hellerman, 1965) to underestimates of the 
computed monthly-averaged wind stress of 10 to 30 percent. 

Recently there has been considerable interest in updating the class- 
ical estimates of Budyko (1963) and Hellcrman (1967). An excellent example 
of this type of work is Bunker's (1976) estimates of the surface heat bal- 
ance over the Atlantic Ocean. The approach used by Bunker differs from 
the classical approach since he computed the turbulent fluxes for each 
simultaneous measurement of surface meteorological variables. The com- 
puted fluxes within a given month are then summed, and finally are divided 
by the total number of observations. Data processing technology has made 
Bunker's approach feasible. However, more than thirty million ship reports 
would have to be processed at the present time to obtain global estimates 
of the turbulent heat and momentum transports. The classical method 
clearly has a number of practical advantages which will be discussed fur- 
ther in our final section. 

The main purpose of this paper is tc provide a basis for determining 
to what extent the classical method of estimating turbulent fluxes of sen- 
sible and latent heat from monthly averaged surface meteorological variables 
can be used for studies of interannual and decadal climatic variability. 

A unique feature of our study is the systematic exploration of the effects 
of boundary layer stability and the physics of the air-sea interface on 
the estimates of the sensible and latent heat fluxes using recent devel- 
opments in boundary layer theory. We will also discuss the importance of 
stability effects for tie estimates of air-sea momentum transfers. Com- 
putations will be based on data taken from ocean weather ships in the North 
Pacific and North Atlantic Oceans since these data are the best available 
long-term records of surface marine meteorological variables. In the final 
section, we will discuss the implications of our results for studies of 
climatic variability. 

2. The Bulk Aerodynamic Calculations 

The well-known bulk aerodynamic formulas for the sensible heat flux 
S and the latent heat flux LE at the air-sea interface are 


and 


5 ‘ » c p c hVVV 


(i) 


LE - pLC E V a (q s -q a ). 


( 2 ) 


where the variables V a , T fl and q fl are the wind speed, temperature, and 
specific humidity of the air at a given height (z_) which is within the 

o 

surface boundary layer. The variables T and q are the sea-surface 
temperature and the saturation specific humidity of air with temperature 

i 

T and sea-level pressure p . In this study we have used z a * 10m and 
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approximated the potential temperature difference (8-8 ) by (T »T ) . 

The variable p is the average air density in the layer from the surface 
to Cp is the specific heat of air at constant pressure; and L 
is the latent heat of evaporation of water. The quantities and 
are the so-called transfer coefficients for sensible heat and evaporation, 
respectively. 

To estimate the fluxes from Eqs, (1) and (2), the transfer coeffi- 
cients must be determined. Although these coefficients may be determined 
in a completely empirical manner using regression analysis, boundary layer 
similarity theory provides a rational, physically motivated way to minimize 
the dependence of and on empiricism. In particular, similarity 
theory allows the inclusion of the effects of the stability of the atmo- 
spheric surface boundary layer; the measure of stability being the Monin- 
Obukov length which is defined as 


L = - 


(t/p) 


3/2 


£ w'T’ 


( 3 ) 


where t is the magnitude of the surface wind stress (see section 7), T 
is the virtual temperature, and g is the gravitational acceleration. Here the 


overbar indicates the Reynolds averaging of the turbulent product w'T^ where 
w is the vertical air velocity. Similarity theory also allows the effects of 
the roughness of the sea -surface to be parameterized by a roughness length z Q < 
We have chosen to evaluate the transfer coefficients using the 
method proposed by Liu et al. (1979) with the exception that we have 

*The variables V , T , T , q and q are all assumed to be averages over a 

3 S 3 S 3 

time-scale at least as long as the lifetime of boundary layer convective 
elements. 
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neglected the difference between the "observed” sea-surface temperature 
and the skin temperature of the water at the air-sea interface. The 
effects of boundary layer stability are included using the Businger-Dyer 
model (see Businger et al., 1971 and Dyer and Hicks, 1970). Also the 
Liu et al. model determines the roughness length z q and the corresponding 
lengths for the profiles of temperature and specific humidity, z T and z ; 
the lengths are not assumed to be equal. There is no free convection 
limit for unstable conditions as considered by Deardorff (1972); however, 
calm conditions with large unstable air-sea temperature and humidity 
differences are rare over the oceans. 

The dependence of and Cg on wind speed and air-sea temperature 
differences computed by the Liu et al. method are shown in Fig. la and lb 
for a relative humidity of 70%. A comparison of these results show that 
for any given wind speed and temperature difference, and Cg differ by 
less than 10%. However, at wind speeds less than 10 m/s both and 
Cg are strongly dependent on both wind speed and air-sea temperature dif- 
ferences. At wind speeds greater than 10 m/s, the stability of the 
boundary layer has very little effect on the transfer coefficients and 
the dependence on wind speed is weak. 

The transfer coefficients proposed by Liu et al. (1979) agree well 
with those of Rondo (1975), and give a reasonable fit of the currently 
available data from the surface boundary layers of the ocean and atmo- 
sphere. However, as we will show in sections 4 and 5, the conclusions of 
this study do not depend crucially on the choice of transfer coefficients. 
This may seem paradoxical to the reader given the importance of the trans- 
fer coefficients in obtaining reliable estimates of the fluxes. However, 
this study is concerned with the methodology of estimating monthly-averaged 
and long-term climatological fluxes from the bulk aerodynamic equations. 

We believe our conclusions concerning the choice of method for obtaining 
flux estimates are valid for any physically reasonable choice of transfer 
coefficients. 

3. Two Estimation Methods for Monthly-averaged Heat Fluxes 

In this section we distinguish between the sampling method and the 
classical method for computing monthly-averaged transfers of sensible and 


latent heat at the sea surface from the bulk aerodynamic equations. In 
the sampling method, fluxes arc estimated from the sample means of S and 
LE in Eqs. (1) and (2) for a given calendar month and geographic area. 

These sample means are defined as 

5 ■ ‘VhVVV- < 4 > 

and 

LE - PL C^V a (q s -q a ), (5) 

where the averaging operation defined by the overbid is 

o - i c )„. m 

and N is the number of simultaneous measurements of T , q , T , and V . 

It should be noted that p, c^, and L are not included in the averaging 
operation since these variations either play a negligible role in the 
surface flux estimates or can be incorporated using linear combinations 
of Sf and LE (e.g.. Brook, 1978). 

The classical method (see Budyko, 1963) uses the following expressions: 


S - ’ c r C H V a <VV ■ 

(7) 

LE ■= H. C E V a (q s -q a ) , 

00 

w 


. where C L1 and C c are obtained by substituting the sample means of T , T , 
He as 

q a , q g , and V 0 into the algorithm for and Cg. 

If measurements of surface meteorological variables are always simul- 
taneous and dense in space and time, and if the resources for data process- 
ing are adequate, the sampling method provides the best estimate of the 
monthly- averaged flux. However, for non- simultaneous and sparse data sets, 
the classical method may have both practical and theoretical advantages 
as discussed in section 8. 

4. Data 

To quantitatively test the differences between the two methods, we 
have analyzed three-hourly marine meteorological observations from nine 
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ocean weather ships. The study was limited to the 25 year period from 
1948-1972 and the ships were selected to give good geographical coverage. 
Table 1 shows the ship positions and periods for which observations were 
available. The position of each ship is also shown in Fig. 2 with the 
ocean surface current regions from Sverdrup et al. (1942). As can be seen, 
the ships have sampled both eastern and western middle-latitude oceanic 
regions including the strong western boundary currents of the Kuroshio and 
Gulf Stream. Six of the ships have nearly continuous records for the en- 
tire 25 -year period. 

r 

We processed the data to eliminate any observations that were obviously 
in error or in which the weather ship was not within 1 degree of the posi- 
tion given in Table 1. We also eliminated data when simultaneous obser- 
vations of wind speed, wind direction, sea level pressure, and air, sea 
and dew point temperature were not available. 

Table 1. Ocean weather ship stations studied. 


Station 

Latitude 

Longitude 

Period Studied 

A 

62°N 

33°W 

1948-1971 

B 

56^°N 

51°W 

1948-1972 

D 

44°N 

41°W 

1950-1972 

H 

36°N 

70°W 

1950-1954 

J 

52^°N 

20°W 

1951-1971 

N 

30°N 

140°W 

1948-1972 

P 

S0°N 

145°W 

1950-1972 

V 

34°N 

164°E 

1956-1971 

X 

39°N 

153°E 

1948-1953 


5. Long-term Climatological Estimates 

In this section we see how closely the classical method approximates 
the sample mean of flux estimates from temporally dense observations taken 
within a given calendar month for a period of ten to twenty years. The 
principal result is shown in Fig. 3 for January and July. The classical 
method agrees with the method of sample means to within a few percent in 
both summer and winter and in all climatological regimes of the subtrop- 
ical and middle-latitude oceans which we have investigated. The only 


exceptions to this conclusion occur when the fluxes arc small. However, 

when the fluxes are small, percentage comparisons are meaningless and 

. -2 

Fig- 3 shows the agreement to be within a few W m in these cases, 

The authors believe the above conclusions are not likely to be 
changed by any physically reasonable choice of transfer coefficients. 

To demonstrate this point, the fluxes shown in Fig. 3 were recomputed 
with the assumption that sensible heat, water vapor, and momentum are 
exchanged between the ocean and the atmosphere in the same manner. (Thus, 
in the Liu et aln (1979) algorithm z T * z^ * z Q .) As shown in Fig. 4, 
the transfer coefficient for sensible heat using this unrealistic assump- 
tion has an extreme dependence on wind speed. The curves for Cg are simi- 
lar. 

The heat fluxes for January with this assumption are shown in Fig. 5. 
Again, the discrepancies between the classical method and the sample mean 
are only a few percent for the latent heat fluxes and about ten percent or 
less for the sensible heat fluxes. This result implies that the difference 
between Budyko' s (1963) flux estimates and Bunker's (1976) estimates for 
the Atlantic Ocean are due primarily to data availability, data reliability 
and the use of different transfer coefficients. The differences are not 
due to the fact that Budyko chose the classical method and Bunker chose the 
sampling method. 

Fig. 6 shows the comparison of the classical and sampling methods using 
only nighttime observations for January. The conclusions based on Fig. 3 
are unchanged. Diurnal effects, whether they are real or instrumental, do 
not lead to significant differences between the two methods. 

We now consider the importance of including boundary layer stability 
effects on long-term climatological estimates of the sensible and latent 
heat fluxes. The heat fluxes with the neutral transfer coefficients from 
Fig. 2 were recomputed using the sampling method. Fig. 7 shows a compar- 
ison of these results with the stability-corrected sailing method results. 
The figure shows differences of approximately ten percent. Budyko (1974) 
reached a similar conclusion by using a bulk Richardson number as the sta- 
bility parameter for the boundary layer. 
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6. Monthly Averaged Estimates 

The results from the previous section suggest that the classi 
method may also be successful in monitoring the interannual variability 
of the sensible and latent heat fluxes, To examine this possibility we 
calculated the heat fluxes for each January and July at all of the ocean 
weather ships using both the classical method and the sampling method. 

Since the results of the comparison are almost independent of the choice 
of month or ship, we have chosen to present only the January results at 
ships D and N representing significantly different heat flux regimes. 

Fig, 8a shows the comparison for Ship D. Although the average heat 
fluxes are large (see Fig. 3), the interannual variability of the monthly- 
averaged fluxes is also large. Despite the observed variability, the 
classical method still gives estimates of the heat fluxes which are within 
ten percent of the values obtained by the sampling method. These results 
are consistent with the conclusions of Fissel et at. (1977) for Ship P. 

They performed cross-spectral analysis on a series of two-year segments 
of data at Ship P with constant transfer coefficients and found that dis- 
turbances with periods of less than a month contributed only ten percent of 
the biennially-averaged sensible heat flux and only five percent of the 
biennially-averaged latent heat flux. 

The results in Fig. 8b for Ship N, show that the interannual varia- 
bility of the heat fluxes in subtropical regions is also well resolved by 
the classical method. The classical estimate of the latent heat flux is 
systematically higher than the sampling estimate, but still within ten 
percent. 

The comparisons in this section were based on the transfer coeffi- 
cients shown in Fig. 1. However, we also computed the monthly-averaged 
fluxes with the unrealistic transfer coefficients shown in Fig. 4. Since 
these coefficients have a stronger dependence on wind speed, the sampling 
and classical' methods are not in such close agreement. However, the 
average disagreement for the 25 year period was still less than 10% and 
never exceeded 20% in any individual month. 

Finally, we wish to point out that the transfer coefficients them- 
selves are subject to uncertainties of approximately 20 to 30% (e.g., Pond 
et at., 1974), In view of this uncertainty, the classical method is a 
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practical and reliable alternative to the method of sample means. The 
uncertainty underscores the need for independent methods for obtaining 
the heat fluxes such as the budget approach of Oort and Vonder Haar (1976). 

7. Climatological and Monthly- Averaged Wind Stress Estimates 

Although the classical method is successful for the estimation of 
sensible and latent heat fluxes, it cannot be used to compute the wind 
stress. As discussed in Hellerman (1965) a wind rose with both direction 
and speed categories is necessary to accurately determine the wind stress. 

In this section we re-examine He Herman's result using the weather ship data 
and also examine the effect of stability on the computation. 

The wind stress vet, *or , x, can be determined from the bulk aerodynamic 
expressions 

T » * " C D V a “x < 9 > 


Here t x and x^ are the east-west components and north-south components of x, 
respectively; u x and are the components of the wind;, and Cp is the drag 
coefficient for momentum, The Liu et al. (1979) algorithm was used to com- 
pute Cp. For a relative humidity of 70%, the dependence of Cp on the wind 
speed is shown in Fig, 9 for several air-sea temperature differences. 

The behavior of Cp is similar to the behavior of and C £ (see Fig. 2) 
except that Cp increases significantly at wind speeds larger than 10 m s" 1 . 

As described in Liu et al. (1979) momentum can be transferred by pressure 
forces acting on the wind as it flows over the waves as well as by frictional 
stresses at the surface. The transfer by pressure forces tends to be strong- 
ly dependent on wind speed. However, molecular diffusion is the only mech- 
anism available for the transfer of sensible and latent heat fluxes at the 
air-sea interface. Thus, the dependence of and on wind speed tends 
to be less than the dependence of Cp for moderate to strong winds. 

Using the Liu et al. algorithm with the stability correction, we com- 
puted the climatological x and x at each weather ship with the sampling 
method and with a wind rose with an average wind speed for each of the 16 
direction categories. Fig. 10 shows that the "direction-only" wind TOse 
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consistently underestimates the wind stress when |t | or j t | is greater 

7 * y 

than 0.5 dynes /cm . This relative error can be as large as 25% (see t 
for Ship B). Thus, Heller-man' s wind rose conclusions for neutral coeffi- 
cient wind stress tlso applies for stability-corrected coefficients. 

Therefore, the wind stress should be computed either by the sampling method 
or by a wind rose with several speed categories for each direction. 

To clarify the importance of stability we recomputed. the wind stress 
with the sampling method for a neutral Cp. The results are shown for the cli- 
matological tg and Ty in Fig. 11. The importance of stability is small and 
affects the results by less than 10% except for small values of t and t . 

A y 

If we reexamine Fig. 7, we see the effect of stability generally produces 
a smaller relative error for wind stress than for sensible and latent heat 
fluxes. 

The effect of stability on t and t was also examined on an inter- 

x y 

annual basis for each ship for January and July. The effect of stability 
was found to be small in each ease. This comparison is shown in Fig. 12 
for the January r x at Ship D. It is apparent that for many applications 
the atmospheric surface boundary layer may be considered to be neutral when 
computing wind stresses on time scales of one month or longer over the open 
ocean. 

The small errors caused by the use of neutral coefficients can be 
further reduced by substituting monthly-mean values of (T s “T a ) and (q g ~q a ) 
into the transfer coefficient algorithms in place of the air-sea temperature 
and humidity differences at each observation time. Fig. 13 shows a compar- 
ison between t x computed with monthly-mean temperatures and humidities and 
t computed using temperatures and humidities at each observation time at 

A 

Ship D during successive Januaries. The errors of less than few percent 
are typical of the results we found at every ship position. 

8. Summary and Discussion 

We have shown that monthly-averaged values of middle-latitude and sub- 
tropical surface marine meteorological data can be used in bulk aerodynamic 
formulas to obtain the monthly-averaged surface fluxes of sensible and 
latent heat. These results confirm and extend the work of Budyko (1963) 
and his coworkers by including recent developments in boundary layer 
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similarity theory to give a more realistic representation of the physical 
processes which are important to air-sea transfers of sensible and latent 
heat. 

Our results also confirm that neutral transfer coefficients give 
reliable wind stress estimates over the ocean. If greater reliability 
is desired, this is easily achieved through the use of monthly-averaged 
air-sea temperature and humidity differences for the determination of the 
drag coefficient. 

Similar results for heat flux and wind stress estimates should be 
obtained in tropical regions. Although this study does not include trop- 
ical marine data, our conclusions are valid for the small Bowen ratios 
(S/LE) and the moderate to low wind speeds which characterize large portions 
of the tropical oceans (see for example Fig. 3b). 

We believe that the classical method of using wind roses and monthly- 
averaged state variables to estimate heat fluxes and wind stresses has a 
number of advantages over computing weighted sample means, These may prove 
to be valuable in describing and monitoring interannual and decadal climate 
fluctuations. 

a- Data processing requirements . The computational effort required 
to process monthly-averaged data from ship-of-opportunity and 
ocean weather ships is much less than the effort needed to process 
all the individual surface reports. Therefore, the revision of 
heat flux estimates with an imporved transfer coefficient or up- 
dated mean field is a relatively simple task with the classical 
method. 

b* Statistical reliability . The classical method may provide a more 
reliable estimate of the heat flux when the number of observations 
is small. Since both the covariance and the correlations of V 
with either (T g -T a ) or (q g -q a ) are known to be relatively small, 
the non-linearities in the flux estimates using the sampling method 
may lead to more statistical uncertainty in estimating the long- 
term average heat fluxes than the classical estimators. 
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c. Utilization of heteorogeneous data . The classical method may 
provide a simple way of combining non-simultaneous data from a 
variety of observing platforms. For example, satellite data does 
not yet provide accurate estimates of the air-soa temperature and 
humidity differences, but it may soon give useful information on 
the surface wind speed. The use of the classical method removes 
the restriction that the observations of wind-speed temperature 
and humidity fields be simultaneous and thus provides a straight- 
forward way to incorporate the satellite observations into the 
heat flux estimates. 
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Figure Legends 

Fig. la. The transfer coefficient for sensible heat from Liu et at. (1979) 
as a function of wind speed at 10m above the sea surface for various air- sea 
temperature differences (AT»T -T ). The neutral surface boundary layer curve 
is also indicated. All curves assume a boundary- layer relative humidity of 
70%. 

Fig. lb. Same as Fig. la except for latent heat. 

♦ 

Fig. 2. Positions of ocean weather ships used in this study superimposed on 
ocean surface currents from Sverdrup et al> (1942) . 

Fig. 3a. Climatological sensible and latent heat flux estimates for January 
for each ship computed by the sampling method (abscissa) and the classical 
method (ordinate). The dashed lines indicate relative errors of 10%. See 
text for details. 

Fig. 3b. Same as Fig. 3a except for July. 

Fig. 4. Transfer coefficients from Liu et at. (1979) for which z Q = z, r = z^ 
(see also Fig. la). 

Fig 5. Climatological sensible and latent heat fluxes as in Fig. 3a but 
computed with transfer coefficients for which z q = = z^. 

Fig. 6. Climatological sensible and latent heat fluxes as in Fig. 3a except 
that only nighttime observations are used. 

Fig. 7. Climatological sensible and latent heat fluxes as in Fig. 3a except 
that values along the ordinate are computed by the sampling method with neutral 
transfer coefficients. 

Fig. 8a. Monthly averaged sensible find latent heat fluxes computed by the 
sampling and classical methods for each January at Ship D. 


- 17 - 


Fig, 8b. Same as Fig. 8a except at Ship N. 

Fig. 9. Same as Fig. la except for the drag coefficient from Liu et at. 
(1979). 

. Fig. 10. Climatological wind stresses computed from the sampling method 
(abscissa) and a "direction-only" wind rose (ordinate). See also Fig. 3a. 

Fig. 11. Same as Fig. 10 except that the ordinate indicates the wind stress 
computed from the sampling method with neutral coefficients. 

Fig. 12. Monthly- averaged t x computed by the sampling method for stability 
corrected and neutral momentum transfer coefficients for earh January at 
Ship D. 

Fig. 13. Same as Fig. 12 except t x (AT, Aq) was computed from stability 
corrected momentum coefficients defined by the monthly-averaged air-sea 
temperature and humidity differences, AT = (T g -T a ) and Aq = 


Transfer Coefficient for Sensible Heat 
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The transfer coefficient for sensible heat from Liu et al. (1979) as a function of 
wind speed at 10m above the sea surface for various air-sea temperature differences 
(&T= T s - T a). The neutral surface boundary layer curve is also indicated. All curves 
assume a boundary- layer relative humidity of 70%. 















Heat Fluxes (W/m 2 ) for January 


Climatological sensible and latent heat flux estimates for January 
for each ship computed by the sampling method (abscissa) and the 
classical method (ordinate). The dashed lines indicate relative errors 
of 10%. See text for details. 
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Fig. 5. Climatological sensible and latent heat fluxes as in Fig. 3a 
but computed with transfer coefficients for which z Q * Zj * 
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Wind Stresses (dyn/cm 2 ) for January 



Fie 10 Climatological wind stresses computed from the sampling method 

(abscissa) and a ''direction-only" wind rose (ordinate). See also 

Fig. 3a. 
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(Neutral) 


Wind Stresses (dyn/cm 2 ) for January 



Fig. 11. Same as Fig. 10 except that the ordinate indicates the wind stress 
computed from the sampling method with neutral coefficients. 






r x for January at Ship D 



Fig. 13. 
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„ «»rpnt t f AT Ao) was computed from stability . 

corrected Lnentum coefficients defined b/Jhe^hlyj^ag^T. 

sea temperature and humidity differences, A T -(T S -T a ) and Aq iq s q a J 


— — «Bwern~— .. 



